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ABSTRACT 

We present results from the Suzaku observation of the microquasar GRS 1915+105 performed during the 
2005 October multi wavelength campaign. The data include both stable state (class x) and limit-cycle oscillation 
(class 9). Correct interstellar absorption as well as effects of dust scattering are fully taken into account in 
the spectral analysis. The energy spectra in the 2-120 keV band in both states are all dominated by strong 
Comptonization of disk photons by an optically thick (r w7-10) and low temperature (T^ «2-3 keV) hybrid 
plasmas containing non-thermal electrons produced with 10-60% of the total power input. Absorption lines 
of highly ionized Fe ions detected during the oscillation indicate that a strong disk wind is developed. The 
ionization stage of the wind correlates with the X-ray flux, supporting the photoionization origin. The iron-K 
emission line shows a strong variability during the oscillation; the reflection is strongest during the dip but 
disappears during the flare. We interpret this as evidence for "self-shielding" that the Comptonizing corona 
becomes geometrically thick in the flare phase, preventing photons from irradiating the outer disk. The low- 
temperature and high luminosity disk emission suggests that the disk structure is similar to that in the very high 
state of canonical black hole binaries. The spectral variability during the oscillation is explained by the change 
of the disk geometry and of the physical parameters of Comptonizing corona, particularly the fractional power 
supplied to the acceleration of non-thermal particles. 

Subject headings: accretion, accretion disks — stars: individual (GRS 1915+105) — techniques: spectroscopic 
— X-rays: stars 



1. INTRODUCTION 

GRS 1915+105 is the brightest microquasar in our Galaxy, 
providing us with a unique opportunity to study the accretion 
flow onto a blac k hole a t high fractions of Eddington ratio 
(for a review, see iFender&Belloni 2004). This source, dis- 
covered in 1992 with the WATCH instrument on GRANAT 
(^astro-Tirado et al. 1992), has been persistently active up to 
now unlike usual soft X-ray transients. GRS 1915+105 was 
recognized to be a super luminal source by VLBA o bserva- 
tions of radio outbursts ( Mirabel & Rodriguez! 1 1994b . mak- 
ing it a target of great importance for studying the formation 
mechanism of relativistic jets. From the kinetics of the ra- 
dio jets, their intrinsic speed and inclination are determined 
to be 0.92c -0. 98c (c is the light spee d ) and 66° -70°, re - 
spectively (Mir abel & RodrfguezI [19941: iFenderet alj|1999l) . 
The near infrared observations have revealed that the sys- 
tem consists of a black hole with a mass of 14 ±4 Mq and 
a K 111 type companion star with an orbital period of 33.5 
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days (i Greiner. Cubv. & McGaughreanll200l1) . This indicates 
a huge size of the accretion disk with a continuously high 
mass transfer rate from the companion, which may account 
for many unique features of this source compared with nor- 
mal black hole binary systems. 

The source has been intensively studied at high energies by 
many observatories on different occasions. RXTE/PCA obser- 
vations revealed that GRS 1915+105 often exhibits dramatic 
temporal and spectral variations, so-called limit-cycle oscilla- 
tions, unique feat ures in accreting stel lar mass black holes in 
our Galaxy (e.g.. iBelloni etanil997al) . iBelloni et alj dlOOOh 
found the instantaneous spectral state of GRS 1915+105 can 
be divided into states. States A, B, and C, where the source 
undergoes frequent transition or stays stable. Based on these 
transition patterns, they classified them into 12 Classes, in- 
cluding 2 "stable" classes (Class x <j>) with little vari- 
ability on a time scale longer than ^ 1 sec. This behavior is 
quite different from normal black holes observed in canonical 
states, the low/hard, intermediat e (or very h igh), and high/soft 
states (iTana ka & Lewin (1995); see also Homan & Bellonil 
(120051) and Remillard & McClintock (2006) for more recent 
classification). The correspondence of States A, B, and C of 
GRS 191 5+105 to these s tates is not completely understood, 
however dReig et alj|2003h . 

Theoretically, the limit-cycle oscillation can be explained 
by thermal instability in the inner part of the accretion disk 
under high mass accretion rate, triggering limit-cycle tran- 
sitions between two stable branches in the surface den- 
sity vs M plane: the standard disk and optically-thick ad- 
vection dominated accre tion flow, so-called the slim disk 
dAbramowicz et al1ll988h . In fact, two-dimensional hydro- 
dynamic simulations are successful to reproduce the limit cy- 
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cle behaviors of GRS 1915+105 qualitatively jOhsuga 120061) . 
However, direct observations that trace the evolution of the 
disk structure both in the oscillation and stable state are still 
insufficient to be compared with model predictions, and their 
interpretations may be subject to large uncertainties since un- 
derstanding of the origin of the X-ray emission is not fully es- 
tablished. Many early works relied on the spectra l modelling 
by a multi-color disk (MCDi' Mrtsuda et al.|[l984l) model plus 
a power law, as applied to other black hole binaries. There is 
no guarantee that this model always holds to GRS 1915+105, 
in particular when the power law component dominates the 
entire flux. 

Later observa t ional works ( Zdziarski et al.l 120011: 
iDone et alj l2004t lUeda et all 120091) have suggested fliat 
the X-ray spectra of GRS 1915+105 at least in some states 
are dominated by an optically-thick Comptonization of 
the disk photons off low energy electrons, by applying 
physically more reali stic models tha n the canonical MCD 
+ power law model. IZdziarski et alj (|2001, 2005) analyzed 
broad band spectra observed with RXTE/PCA, HEXTE, 
and CGRO/OSSE, and found that they can be described by 
Comptonization from thermal/non-thermal hybrid plasmas. 
Such low temperature, optically-thick Comptonization may 
be common characteristics in accretion flows at high Edding- 
ton fraction s, which may apply to those onto supermassive 
black holes jMiddleton et alj|2009l) . 

Observations of local spectral features such as iron-K emis- 
sion line and absorption lines give key information to reveal 
the structure of accretion disks. From high resolution spec- 
tra (CCD or the HETGS instrument on Chandra) of GRS 
1915+105, absorption lines from highly ionized ions have 
been clearly det ected jKotani et al.1 120001: iLee et alj 120021: 
lUeda et aLll2009l) . This indicates the presence of strong disk 
wind most probably occurring at r '--^ lO^rg (rg = GM/c^ is 
the gravitational radius where G, M, and c being the gravi- 
tational constant, mass of the black hole, and light velocity, 
respectively), which carries a huge amount of accreting ga s 
outward the system (lUeda et al.ll2009l: iNeilsen & Le3 120091) . 
On the other hand, the emission line profile gives constraints 
on the geometry between the continuum emitting region and 
reflector, most probably the accretion disk. An iron-K emis- 
sion line has been detected from GRS 1915 +105, whose shape 
and intensity seem to depend on the states jKotani et al.ll200(A 
INeilsen & Leell2009l) . Since its equivalent width is not large 
(< 50 eV) in GRS 1915+105, we need both good spectral 
resolution and large effective area to accurately measure the 
profile, in particular to trace the change during oscillations. 

Suzaku (Mitsuda et al. 2007), the 5th Japanese X-ray satel- 
Ute, observed GRS 1915+105 on 2005 October 16-18 (UT 
throughout the paper) as a part of the science working group's 
observation program. The large effective area with good en- 
ergy resolution in the 0.2-12 keV band, together with the si- 
multaneous coverage of the 10-600 keV band with unprece- 
dented sensitivities, give us the best opportunities to reveal 
the origins of the X-ray emission and disk structure of GRS 
1915+105 both in stable and oscillation states. In our work, 
we refer to the latest results on the metal abundances of inter- 
stellar (plus circumstellar) gas toward GRS 191 5+105 as de- 
termined by Chandra/HETGS (lUeda et al.ll2009h . In the anal- 
ysis, we also take into account the effects by dust-scattering, 
which have been ignored in most of previous works. We show 
that these are important for accurate modelling of the spectra. 
§ 2 describes the observation and data reduction, and § 3 the 



light curves and states of GRS 1915+105 in our observations. 
The spectral analysis is presented in § 4. We discuss the in- 
terpretation of our results in § 5. The conclusions are summa- 
rized in § 6. In our paper, we assume the distance of D = 12.5 
kpc and inclination angle of / = 70° unless otherwise stated. 
Around the epoch of our Suzaku observations, a large mul- 
tiwavelength campaign was conducted involving other space 
and ground observatories, whose results will be reported in 
a separate paper (Ueda et al., in prepara tion). For the pre - 
liminary results of the campaign, refer to Ueda et al. (2006^ 
The INT EGRAL light curves are published in .Rodriguez et aP 
(l2008ah . 

2. OBSERVATIONS AND DATA REDUCTION 

Suzaku observed GRS 1915+105 from 2005 Oct 16 16:42 
to Oct 18 23:16 for a net exposure of «80ksec. Suzaku carries 
four sets of X-ray Telescope (XRT) each with a focal plane X- 
ray CCD camera, the X-ray Imaging Spectrometer (Front-side 
Illuminated XlSs, XlS-0, XlS-2, and XlS-3: Back-side Dlu- 
minated XIS, XlS-1), and a non-imaging collimated instru- 
ment called the Hard X-ray Detector (HXD), which consists 
of the silicon p-intrinsic-n photo-diodes (PIN) and scintilla- 
tion counters made of gadolinium silicate crystals (Ce-doped 
GdiSiOs : GSO). The XIS, PIN, and GSO simultaneously cov- 
ers the energy band of 0.2-12 keV, 10-60 keV, and 40-600 
keV, respectively. The target was observed at the so-called 
XlS-nominal position. 

The XlSs were operated in the normal clock mode with the 
2x2 editing mode for XIS -0, 2, and 3, and 3 x 3 for XlS-1. To 
minimize pile-up (double photon events) due to the brightness 
of this source, the 1/8 window option was adopted for XIS- 
0, XlS-2, and XlS-3, which reduced the exposure time per 
frame from 8 sec to 1 sec. XlS-1 was operated with the 0.1 
sec burst plus 1/8 window option, although the 1/8 window 
option was turned off for XlS-1 before Oct 18 00:26, due to 
an accidental error in the satellite operation. We utilize XIS 
data taken with different modes according to the purpose of an 
analysis. The XlS-1 events are used to derive (1) unsaturated 
light curves with a resolution of 1 sec (after Oct 18 00:26) 
and 8 sec (before then) and (2) continuum spectra with little 
effects from pile-up, while we use XlS-0, 2 and 3 data to make 
the spectra with high photon statistics in the 5-10 keV band 
for detailed study of the iron-K features. 

We analyze the data using the HEAsoft version 6.4.1 pack- 
age, utilizing the products version 2.0.6.13 processed by the 
pipeline processing team and calibration database (CALDB) 
released on 2008 Dec. 3. For the analysis of the XIS data, we 
exclude events suffering from telemetry saturation, based on 
the GTl filters provided by the XIS teanfl. We do not subtract 
the non X-ray background nor the cosmic X-ray background, 
which are negligible compared with the bright source counts. 
For the spectral analysis of the XlS-1 data taken with the full 
window mode before Oct 18 00:26, we subtract out-of-time 
events as background by utilizing events in an outer region in 
the same CCD chip. In the extraction of the spectra of XlS-0, 
2, and 3 (i.e., without burst options), we exclude a circular 
region with a radius of 120 pixels around the target position 
to reduce the effects of pile-up events. It is known that the 
attitude of the Suzaku satellite exhibits an orbital variation by 
^0.5 arcmin, causing periodic shifts of the target position in 
the focal plane. Hence, we check the central position of the 
source and define an event extraction region every 200 sec- 

' http://www-cr.scphys.kyoto-u.ac.jp/menber/hiroya/xis/gtifile.html 
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Fig. 1.— The light curves of GRS 1915+105 in the 1-10 keV and 10- 
60 keV band obtained with XIS-1 and HXD/PIN, respectively. The data are 
binned in intervals of 16 seconds. The red aiTows denote the time regions 
used for the spectral analysis of the stable and oscillation states. 

onds. 

For the HXD, we use the best modelled background files 
provided by the HXD tearrFI. The reproductivity of the NXB 
spectra with a 10 ksec exposure is known to be 1 .4% and 0.6% 
for the PIN (15-40 key ) and GSO (50-100 keV), respectively 
jFukazawa et al.l 120091) . Then, an estimated CXB spectrum 
is added to the NXB spectrum of the PIN with the angular 
transmission function for a uniformly extended emission. The 
CXB is neglected in analysis of the GSO data. 

3. LIGHT CURVES AND X-RAY STATES 

Figure 1 shows the light curves of GRS 1915+105 obtained 
with XIS-1 in the 1-10 keV band and with the HXD/PIN in 
the 10-60 keV during the whole observation epoch. All pho- 
tons observed in the XIS-1 chip are utilized here, but the XIS- 
1 count rates before Oct 18 00:26 are artificially reduced by 
a factor of 0.51 to correct for the difference of the integra- 
tion area from the data taken with the 1/8 window option. For 
convenience, hereafter we define t as the time since 2005 Oc- 
tober 16.0 in units of hours. It is seen that at f ^40 GRS 
1915-1-105 was in an X-ray stable state, except for a possi- 
ble flare-like event seen at f ~ 21.5. From f ~ 40, the X-ray 
flux increased rapidly and entered into "oscillation" that lasted 
till f ~ 56. From t ~ 66, GRS 1915H-105 exhibits a similar 
X-ray flare, followed by oscillation whose variability pattern 
is different from the pr evious one (see Ueda et al...2006. and 
iRodriguez et"al]l2008al for X -ray light curves observed with 
other satellites). Throughout the observation, the source was 
very bright in soft X rays, ^ (1 -3) x 10"^ erg cm"^ s"' in the 
2-10 keV band. 

For comparison with previous studies of GRS 1915-1-105 
we first identify the "Class" and "State" of GRS 1915H-105 in 
our observations, according t o the model- independent clas- 
sification scheme by Belloni 'et al.l ( 120001) . As seen from 
Figure 1, GRS 1915-1-105 continuously stayed in the "stable 
state" in the period of f = 26 - 38, which is similar to the 
plateau state dFoster et al.ll 19961) . before entering into the os- 
cillation phase. Timing and spectral analyses of the data taken 
with the proportional counter arrays (PCA) on RXTE at the 
epoch of f = 23 - 26 reveal that our "stable state" corresponds 

http://www.astro.isas.ac.jp/suzaku/analysis/hxd/pinnxb/tuned/ (PIN) 
and http://www.astro.isas.ac.jp/suzaku/analysis/hxd/gsonxb/ (GSO) 
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Fig. 2.— A blow up of XIS-1 (1-10 keV) and PIN (10-60 keV) light curves 
att = 51.6-52.3 and their hardness ratio with 4 sec bin. The two horizontal 
lines indicate the limits between the oscillation low/med and med/high states, 
respectively. 

to Class X (stable State C) having an extremely soft spectrum, 
acc ompanied with a quasi periodic oscillation (QPO) at fti6 
Hz (U eda et al.ll2006h . 

Figure 2 shows representative X-ray light curves (f = 5 1 .6 - 
52.3) in the oscillation phase taken with the Suzaku XIS-1 
and PIN in the 1-10 keV and 10-60 keV band, respectively, 
together with their hardness ratio. We iden t ify the X-ray vari- 
ability pattern as Class 9 in 'B elloni et al.l (120001) . where the 
transition occurs between State C and State A (soft dip) on a 
time scale of ^10 minutes. It is seen that during oscillations 
the spectral colors significantly softens in the dip phase, as ex- 
pected according to previous studies with RXTE (iBelloni et al.l 
i2000i) . 

4. SPECTRAL ANALYSIS 
4.1. Overview 

For spectral analysis, we define four different states (a) 
"stable state" (Stable), (b) "oscillation high state" (Osc-H), 
(c)"oscillation medium state" (Osc-M), and (d)"oscillation 
low state" (Osc-L), according to time variability pattern and 
intensity level. The stable and oscillation-high stat es con- 
sist of State C data defined bv IBelloni et al.ll2000l) . while 
oscillation-med and low corresponds to State A. The time 
region is taken between t = 26.03-37.88 for (a), when the 
source flux was almost constant (Class x), and t =45.18- 
53.97 for the rest (Class 9). The oscillation high, med, low 
states are defined as when the 10-60 keV PIN count rate is 
above 24 c s"^ between 12-24 c s"', and between 1-12 c 
s"' , respectively, within the above time region. The count rate 
thresholds for the oscillation states are illustrated by two lines 
in the middle panel of Figure 2. Table 1 summarizes the defi- 
nition of the states. We make final good time intervals (GTIs) 
by combining the XIS, PIN, and GSO, to cover exactly simul- 
taneous time regions in each state. The net exposures are 19.8 
ks, 7.2 ks, 0.93 ks, and 2.6 ks for the stable, oscillation-high, 
med, and low state, respectively. 

The spectral analysis we perform consists of two steps, (1) 
analysis of local features in the 5-10 keV band, which deter- 
mines the parameters of a reflection component (see below), 
and (2) that of continuum emission over the 2-120 (or 2-50) 
keV band. As mentioned above, different XIS sensors are 
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TABLE 1 
State Definition 



State 


Start 


End 


Criteria 


Exposure 




(UT) 


(UT) 


PIN (10-60 keV) 


(ksec) 


Stable 


2005 Oct. 17 02:02 


Oct. 17 13:53 




19.8 


Osc-H 


2005 Oct. 17 21:11 


Oct. 18 05:58 


>24 c/s 


7.2 


Osc-M 


2005 Oct. 17 21:11 


Oct. 18 05:58 


12-24 c/s 


0.93 


Osc-L 


2005 Oct. 17 21:11 


Oct. 18 05:58 


1-12 c/s 


2.6 



used for each purpose. For the former step, we use XlS-0, 2, 
and 3 (Fl-XlSs), which were operated without the burst op- 
tions, by excluding the central core of the point spread func- 
tions; this makes us achieve good photon statistics by avoid- 
ing extreme pile-up effects. The spectra of the three sensors 
are summed together with the energy responses. For the lat- 
ter step, we make a simultaneous spectral fit to the spectra 
of XlS-1 (Bl-XlS), PIN, and GSO, covering the 2-120 keV 
band, for the stable and oscillation-high states, while we only 
use the XIS and PIN spectra in the 2-50 keV band for the 
oscillation-med and low states, since the GSO data in these 
states have limited statistics due to the low flux and short 
exposure. The data of XlS-1, operated with the 0.1 s burst 
option, is almost free from pile-up effects and we can uti- 
lize all photons in the point spread function. Thanks to this 
fact, we can best estimate the continuum below 10 keV with 
least uncertainties in the calibration of the energy response of 
the XIS+XRT system, although the photon statistics is poorer 
than the data of the three Fl-XISs. The observed spectra of the 
XlS-1 and HXD in the stable and oscillation state are plotted 
in Figure 3 (a) and (b), respectively. 

Since the determination of the local features affects that of 
the continuum and vice versa, the above two steps are re- 
peated in an iterative way. Firstly, we roughly determine the 
continuum shape from the BI-XIS and HXD spectra. Sec- 
ondly, using the Fl-XISs spectra in the 5-10 keV band, we de- 
termine local iron-K features, including absorption lines and 
an iron-K emission line, which is modelled by a feature from 
a reflection component. In this stage, we fix the parameters 
of the continuum at the best-fit obtained above except for its 
normalization. To take into account a remaining component 
of pile-up events, we add a power law component having a 
positive photon index in the model, whose parameters are set 
free. Thirdly, we again perform the continuum fit of the Bl- 
XIS and HXD spectra, but including the reflection component 
and absorption line features, whose parameters are fixed at the 
best-fit values as determined in the previous fit of the Fl-XISs 
spectra. Detailed description of the continuum model and lo- 
cal features is given in the following subsections. 

4.2. Absorption and Dust Scattering 

In our study, special care is paid on two important effects 
that could affect the modelling of the continuum, (1) inter- 
stellar absorption and (2) dust scattering. It is well known 
that the elemental abundances in interstellar medium in our 
Galaxy are different by location, and are not "Solar", as often 
assumed in spectral analyses of Galactic sources. To model 
the total photo-electric absor ption t oward GRS 1915+105, we 
refer to the results by U eda et al.l ('2009) obtained from the 
C/zanc/ra/HETGS spectrum, in terms of "equivalent" hydro- 
ge n column densities (conve rted from the Solar abundances 
bv lAnders & Grevesse|[T989l) of 2.78 x lO^^ cm'^ for H, He, 
C, N, and O, 6.3 x 10^^ cm'^ for Ne, Na, Mg, and Al, 



5.8 X 10^2 cm-2 for Si, 7.6 x 10^^ cm'^ for S, CI, Ar, and 
Ca, and 10 X 10^^ c m"^ for Cr, Fe Co, and Ni. The absorp- 
tion cross section bv IWilms et al.l ( 120001) is adopted, which is 
available as TBvarabs model in XSPEC. 

We also take into account the effects by dust scattering, 
which are more significant in soft X-rays particularly below 

3 keV. A part of direct emission from a heavy absorbed 
source is scattered out by interstellar dust in the line-of-sight, 
while photons emitted with slightly different angles are scat- 
tered in, thus making dust-scattering halo around the target. 
At first order approximation the scattered-in and scattered-out 
lights cancel each other, although there is time delay before 
scattered-in components reach an observer following the di- 
rect component. Hence, we do need to correct for this effect 
as far as the source spectrum is constant and one can integrate 
all the emission including both direct and halo components. 

This is not the case for GRS 1915+105 observed with 
Suzaku, however, because (1) the halo components cannot 
be fully collected in a limited field of view of the XIS and 
(2) the source is highly variable in the oscillation states. As- 
suming th e same image pr ofile of the dust scattered halo as 
GX 13+1 ( Smit h et al.ll200 2). which has a similar absorption 
column density to GRS 1915+105, we find that the effective 
area of the XlS-1 (full window) for the halo component is 
«60% of that for the direct (i.e., point source) one, and thus 
the scattering effect is not canceled out even when the spec- 
trum is constant. When it is variable, the situation becomes 
more complex; the total spectrum observed in a given time in- 
terval is contaminated by halo components emitted in earlier 
epochs due to their time delay. The time scale of the delay is 
estimated to be ^ fP-Djc, where Q and D is the halo size and 
distance to the source, respectively. Taking 9^2' according 
to the halo profile of GX 13+1, we estimate the delay to be 
^100 hours, much longer than the time scale of the oscilla- 
tions. This means that we approximately observe the spec- 
trum of variable direct emission plus constant dust-scattered 
component in the oscillation states . According to the light 
curve of all sky monitor on RXTE ( Ueda et ani2006l) . the av- 
eraged flux level was almost constant for ^-^15 days before our 
observation epoch, and thus we can reasonably assume that 
the shape of the direct spectrum producing the observed halo 
component is the same as that of the stable state. 

To make the best estimate of the continuum spectra includ- 
ing low energies below ^3 keV, we derive the spectrum of the 
halo component in the stable state, by taking into account the 
fact that 60% of the scattered-in component is contained in the 
observed spectrum. Here we refer to the cross section of dust 
scattering given by lDrain3 (120031) . which is fine- tuned to re- 
produ ce the observed halo intensity of GX 13+1 (ISmith etal J 
120021) at its hydrogen column density dUeda et al] 2005h . For 
the oscillation states, we subtract this component from the ob- 
served spectra and fit them by considering the opacity of both 
dust-scattering and photo-electric absorption. We assume that 
the halo component is the same as in the stable state and is 
constant during the oscillation. The dust-scattered component 
is plotted in Figure 3 (a) and (b). Its fraction in the direct 
(unscattered) component is found to be 6%, 7%, and 15% 
at 3 keV in the oscillation-high, med and low state, respec- 
tively, which increases toward lower energies approximately 
as cx E~^. 



4.3. Iron-K Features 
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Fig. 3. — (a) {left) The observed spectra of XIS-1 (full window mode) and the HXD in the stable state, folded by the detector response, (b) (right) Those of 
XIS-1 (1/8 window mode) and the HXD in the oscillation-high (upper, black including the GSO), oscillation-med (middle, red), and oscillation-low (lower, blue) 
states. The lower curves in the XIS-1 band (green) indicate the estimated contribution of the dust scattered component. 



The spectra of the FI-XISs with rich photon statistics re- 
veal that the iron-K features are different between the four 
states. Figure 4 shows the summed FI-XISs spectra in the 5- 
10 keV band from which the estimated pile-up components 
(see above) are subtracted. The overlaid model is the best- 
fit continuum model with its reflection component, described 
below, from which both iron-K emission and absorption lines 
are excluded to show their significance; the ratio of the data 
to the model is plotted in the lower panel. In addition to a 
deep absorption-edge at 7.11 keV due to the overabundance 
of iron in the interstellar or circumstellar medium toward GRS 
1915-1-105, we detect an emission line centered at 6.5-6.6 keV 
in both stable and oscillation low states, whereas it is very 
weak in the oscillation high and medium states. We also find 
a strong absorption line feature from He-like and H-like iron 
ions at ~ 6.6 and ~ 7.0 keV in the oscillation low state, which 
is weaker in the other states; in particular, no significant ab- 
sorption lines are detected in the stable state. 

We model the iron-K emission line by utilizing the pexriv 
reflection code ( Magdziarz & Zdziarski 1995), which cal- 
culates both reflection continuum and fluorescence line. 
The line profile is further blurred with the diskline kernel 
( lFabianetal.1119891) . by assuming that it originates from the 
accretion disk between the radii Hn and rout with an emis- 
sivity law of . The free parameters are (1) the solid an- 
gle of the reflector Vl, which basically determines the equiva- 
lent width (EW) of the emission line, (2) ionization parameter 
^ = Ljnp- (Tarter et al. 1969) (where L, n, and r is the lumi- 
nosity, density of the reflector, and distance from the emitter, 
respectively), and (3) innermost radius ri,,. We fix an iron 
abundance at the same value as in the interstellar absorption, 
an inclination of 70°, /3 = -3, and rout = 10'' /"g. Note that the 
pexriv model does not properly calculate the additional broad- 
ening by Comptoniza tion in the emission line profile (see 
iBallantyne et al.lf2001l) . Thus, we may slightly underestimate 
the true equivalent width, and hence the reflection strength in 
the stable state, where the ionization parameter is found to be 
log ^ = 2-3. In the oscillation states, the effects are expected 
to be negligible due to the low ionization parameter (log ^ ^ 
2). The absorption lines are modelled by Gaussians with a Ict 
width fixed at 12 eV. The best-fit parameters are summarized 
in Table 2. 



4.4. Continuum Model 

The overall shape of the continuum in the four states is char- 
acterized by a steep power law with a photon index of « 3 at 
^5-50 keV. To find an appropriate model to describe the con- 
tinuum, we begin with fitting the XIS -i- HXD spectra in the 
2-120keV in the stable state, by taking into account the inter- 
stellar absorption and dust-scattering described in § 4.2, and 
the reflection component whose parameters are determined in 
the previous subsection. In addition, to include the opacity 
of the highly ionized gas responsible for the iron-K absorp- 
tion lines, we introduce absorption edges at ss8.6 keV and 
Ri9.0 keV corresponding to those by Fe XXV and Fe XXVI 
ions, respectively. We fix the optical depth according to the 
equivalent width of the correspo nding absor ption line (see Ta- 
ble 2), based on the results by lUeda et al.l (12009) . We first 
apply a canonical model consisting of a MCD component 
and a power law, which turns out to not acceptable (x^/dof 
= 559/354); moreover, the best-fit results are quite unphysi- 
cal because (1) the innermost temperature of the MCD model 
is high (^3 keV) with an extreme small radius (^4 km) and 
(2) the power law component dominates the flux over the en- 
tire band. This confirms that the spectral state is very differ- 
ent from the canonical high/soft stat e of black hole bi naries. 
Similar fitting results are reported by lMuno et al.l ([1999) from 
RXTEI¥CK data of GRS 1915H-105 when the 0.5-lOHz QPOs 
are observed. 

Next, we apply a thermal Comptonization model by 
IZycki. Done. & SmithI (^1999^ with seed photons fr om a MCD 
component, as ado pted by iDone et al.l ( 120041) . lUeda et al] 
(2009), and Vierdayanti et al.l (|2009|). By considering that a 
part of the MCD photons can escape without being Comp- 
tonized, the model is expressed as thComp + diskbb in the 
XSPEC terminology. The innermost temperature of seed pho- 
tons input to the thComp model and that of the MCD com- 
ponent is tied each other We find that, while it gives a 
reasonable fit to the data below ^50 keV, a significant ex- 
cess remains above this energy, indicating the p resence of 
a non -thermal tail as d iscovered with OS SE ( Zdz iarski et al.l 
120011) and INTEGRAL (iRodri guez et al.l2008b ). Accordingly, 
by adding a power law component (and its reflection) to the 
above model, we obtain an acceptable fit (x^/dof = 344/350) 
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Fig. 4. — The FI-XISs spectra in the 5-10 keV band. The soHd curves (red) represent the best-fit continuum model from which an iron-K emission and 
absorption lines are excluded. The lower panel plots the residuals in units of ratio between the data and above model. From left to right and top to bottom, (a) 
stable, (b) oscillation-high, (c) oscillation-med, and (d) oscillation-low state. 



from the broad band spectra in the stable state. This result in- 
dicates that the spectrum can basically be represented with 
Comptonization of disk photons both by thermal and non- 
thermal electrons. The power-law description of non-thermal 
Comptonization is unphysical, however, since it must breaks 
at energies of seed photons in reality. 

As a more physically self-consistent model of Comptoniza- 
tion, we finally adopt the eqpair model developed by Coppil 
(1^999), ins tead of thComp. It was adopted by Zdziars ki et alj 
OOOl.) and lZdziarskietalJ (l2005h to fit the KXTE and OSSE 
spectra of GRS 1915H-105 in Classes 7' ^nd w successfully. 
This model computes Comptonization from a hybrid plasma 
consisting of thermal and non-thermal electrons, where the 
physical processes are self-consistently solved with an en- 



ergy input of accelerated electrons in a background thermal 
plasma. The model parameters are (1) the optical depth for 
scattering r, (2) compactness parameters, defined as Zh.s = 
L'a.^cri/iRm^c^), where Lh.s is a power supplied to the plasma 
and that in seed photons, aj Thomson cross section, R the 
plasma size, nje the electron mass, and (3) its ratio between 
non-thermal and thermal electrons, Int/h, where 1^ = Ini + h- 
The injected non-thermal electrons have a power law energy 
distribution between the Lorentz factors 7inin and 7[nax with 
an index of Finj. In the spectral fitting, we fix = 1000, 
7min = 1-3, 7max = 1000. In this model, the energy balance 
between Compton and Coulomb interactions determines the 
temperature of thermal electrons, which is automatically cal- 
culated from the given parameters. 
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TABLE 2 
Spectral Parameters 



Parameters 



Stable 



Osc-H 



Osc-M 



Osc-L 



Column Densities" (A'h) 



H (10^2 cm-2) 


2.78 


Mg (10^2 cm-2) 


6.3 


Si (10^2 cm-2) 


5.8 


S (1022 


7.6 


Fe(1022 cm-2) 


10 



Continuum'' 



r,„ (keV) 
(n (keV)) 

/h/4 

Int/k 



- inj 



0.53±0.04 

2 2+"" 

Qd+^'^t 
"■'^ -0.02 

0.55±0.02 
8.3±0.3 
3.7±0.2 



'(kmf 
(km)'' 



3001: 
420+^ 



0.79:«;»2 

2.8:«;? 
0.84^62 

n :!Q+0.04 
'-"-0.1 

3.1±0.3 
120+'" 



0.57±0.16 



1.7 



+0.4 



n 4Q+0.14 



a+o.-3 



0.69+; 



210+ 



)+0.10 
-0.07 

0.34±0.09 

10-1,7 
3.1 (fixed) 
150!^" 
360+^^^ 



0.75+ 



+0.16 
-0.07 
0.17±0.05 

10-1.1 
3.1 (fixed) 

190+1* 
390+380 
-'^"-160 



Reflection"^ 



Q/(27r) 


0.11 ±0.02 


< 0.02 


0.10±0.07 


54+0 06 
'^■•'^-0.07 




530+'* 


120 (fixed) 




39!^o 




1000_200 


290(fixed) 


290(fixed) 


200+170 
^^"-140 


(E.W.)' (eV) 


33 


<4 


13 


42 



Iron-K Absorption Features 



£ce„ (keV) 
E.W. (-eV) 



6.76±0.33 

26+7 



£ce„ (keV) 
E.W. (-eV) 



» (fixed) 
<0.9 



6.99+ 



a+il.U'l 
-0.05 
10 1+2'' 



7.01±0.10 7.02±0.02 
15+^ 39+5 



edge 

edge 
^9.0 








0.051 




0.072 


0.11 
0.20 


Luminosity and Flux 


^aoT-ioo^lO'^'^^g^-') 

CT-.oo^(l«''-g^"') 
F2-10'' (10"** ergcm-2 s^') 

fio-so" (10-** ergcm-2s-') 


9.2 
10.4 
1.69 
0.63 


7.6 
15.2 
2.88 
1.03 


3.0 
14.4 
1.92 
0.39 


2.8 
8.8 
1.12 
0.17 


Fitting Information 


Energy Range (keV) 
/ d.o.f 


2-120 
319/351 


2-120 
890/930 


2-50 
114/101 


2-50 
127/128 



Note. — The errors are 90% confidence level for a single parameter. 
" Equivalent hydrogen column densities as determined by Ueda et al] )2009l) , 
in units of Solar abundances I Anders & Grevesse 1989) between the element 
and hydrogen. The abundance ratios within each group of H-He-C-N-O, Ne- 
Na-Mg-Al, S-Cl-Ar-Ca, and Cr-Fe-Co-Ni are fixed at the Solar values. 

Fit with eqpair with a MCD spectrum as seed photons (see text). The com- 
pactness parameter for the seed photons is assumed to be /, = 1000. Non- 
thermal electrons are considered between 7niin = 1-3 and 7max = 1000 with a 
power law index of Finj . 

The normalization of the directly observed MCD component in terms of 
an innermost radius. The distance D = 12.5 kpc and inchnation i = 70° are 
assumed without corrections for the color and boundary condition. 

The innermost radius estimated from both Compton-scattered and unscat- 
tered MCD photons (see text). 

Compton reflection based on the code by iMagdziarz & Zdziarski|[T993) . 
blurred by the diskline kernel with an emissivity law of between the inner 
radius ri„ and outer radius Tout, where /3 = — 3 rout = lO^rg. The accompanied 
iron-K emission is included. 

' The equivalent width of the iron-Ka emission line with respect to the total 
continuum. 

° Intrinsic luminosity of the direct continuum (without the reflection compo- 
nent) in the 0. 1-100 keV band, corrected for absorption. A spherical geometry 
is assumed for the eqpair component, while a disk geometry is for the MCD 
component. 

Observed flux in the 2-10 keV or 10-50 keV band. 



We find the eqpair + diskbb model, with the reflection of 
the former component, provides an excellent description of 
the continuum in all four states. The direct MCD compo- 
nent is required in the stable and oscillation-high states with 
the F-test probability of 7 x lO"-' and 3 x 10"^, respectively. 
Although it is not significant in the oscillation-med and low 
states, we include this component for consistency. Figure 5 
plot the unfolded spectra in units of EI{E), where HE) is the 
energy flux. The Comptonized component, MCD model, and 
the reflection component including the iron-K emission line, 
are separately plotted. The best-fit continuum parameters are 
summarized in Table 2, together with the observed 2-10 keV 
and 10-50 keV fluxes and absorption-corrected 0.01-100 keV 
luminosities for the Comptonization and MCD component, 
where a spherical and disk geometry is assumed, respectively. 
For the oscillation med and low states, we fix Finj = 3.1, the 
best-fit value obtained in the oscillation high state, since it is 
difficult to constrain the slope without the coverage above 50 
keV in these states. 

5. DISCUSSION 

5.1. Results Summary 

With Suzaku, we obtained high quality broad band X-ray 
spectra of GRS 1915H-105 over the 2-120 keV band both in 
the stable (Class x) ™d oscillation states (Class 9), covered 
with CCD energy resolution (FWHM -120 eV at 6 keV) be- 
low 10 keV. The best estimated interstellar absorption and 
dust scattering effects are taken into account in the spectral 
fitting, which have been ignored in most of previous studies of 
GRS 1915H-105. In particular, we caution that without correct 
modeling of the absorption with over iron abundances, one 
can easily reach wrong conclusions on the iron-K emission 
line profile that requires very careful analysis of the contin- 
uum shape. The best-fit models in units of EI(E) corrected for 
the interstellar absorption and scattering are plotted in Fig- 
ure 6 for the four states. Below, we interpret our results based 
on the iron-K features and continuum model. We find that the 
iron-K features in the oscillation low state are very similar to 
those in the "soft state" (stead y State A) of GRS 1915H-105 
detected by lUeda et all (120091) using ChandralWEIGS. This 
may be expected, since our oscillation low state mostly corre- 
sponds to State A, though not in a steady state. 

5.2. Change of Reflection 

We detect significant iron-K emission lines, direct evi- 
dence for the reflection, in the stable, oscillation-med, and 
oscillation-low states. The line profiles indicate that the re- 
flector is moderately ionized (log ^ — 1 - 3) and is located far 
away from the black hole as constrained by its "narrow" line 
width; it is most probably the outer parts of the accretion disk, 
whose inner radius is estimated to be ^ 300rg by assuming an 
emissivity law of r~^ . This does not mean that the disk must 
be truncated at that location. If the scale height of the central 
source is smaller compared with that of the outer disk, then the 
emissivity law is better expressed with an index lower than 3, 
which would yield a much smaller inner radius. The absence 
of a relativistic broad iron-K line is consistent with the picture 
that the innermost part of the disk is covered by an optically 
thick hot plasma, smearing out any emission features due to 
strong Comptonization. 

During the limit-cycle oscillation, the EW of the iron-K 
emission line (see Table 2) is the largest in the oscillation low 
state. It decreases with X-ray flux, and the line disappears 




Fig. 5. — The unfolded spectra in units of EI(E) obtained from tlie XIS-1 and HXD data (red. crosses) plotted with seperate contribution of the model; from 
upper to lower, the total (solid, green), Comptonized component (dot, blue), direct MCD component (dash, cyan), and reflection component (dot-dash, black). 
From left to right and top to bottom, (a) stable, (b) oscillation-high, (c) oscillation-med, and (d) oscillation-low state. 




Energy (KeV) 



Fig. 6. — The best-fit models in units of EI(E) corrected for the interstellar 
absoiption and dust-scattering in the four states. From upper to lower (at 
20 keV), oscillation-high (red), stable (cyan), oscillation-med (green), and 
oscillation-low state (blue). 

in the oscillation high state. The change of the ionization pa- 
rameter, determined by the X-ray luminosity, is expected to 
be within a factor of 3 during the oscillation. Hence, we must 



attribute this change mainly to the variability in the solid an- 
gle of the reflector as seen by the central source. In the stable 
state, the solid angle is small fl/lTT « 0.1 as well, similar to 
that found in the oscillation med state. 

We interpret that this is the first observational evidence 
for "self-shielding" effects of the inner part of the accretion 
disk. The dynamical structure of outer parts of the accre- 
tion disk, r ^ 300rg, cannot be changed on a time scale of 
the hmit-cycle, ~1000 sec, mu ch shorter than that of viscos- 
ity dShakura & Sunyaevl ll973^. Therefore, the only possible 
mechanism that can dramatically reduce the observed reflec- 
tion strength in the oscillation-high state must be related to the 
geometry in the central source; most of the continuum X-ray 
emission seen by an observer (i.e., with / = 70°) does not reach 
the outer parts of the disk in this state. As we describe later, 
the inner accretion disk is covered by hot corona responsible 
for Comptonization. Our observations imply that the inner 
part of the disk becomes geometrically thick in the oscillation 
high state, and Comptonized photons are shielded by the sur- 
rounding cool region of the expanded disk when viewed at a 
very high inclination angle from the outer disk. The evolution 
of the disk geometry is in acco rdance with the oretical predic- 
tion of the slim disk ( Abramo wicz et al.ll 1988b and numerical 
simulation (e.g.. lOhsuga Il2006i) . 
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5.3. Disk Wind 

As have been revealed by pre vious studies ([Kotani et al] 
IMKllLee et al.l2002tllJedaer ar2009:'Neil sen & LeeP 20m 
absorption lines of highly ionized ions indicate the presence 
of a disk wind in the line-of-sight between the continuum 
emitter and the observer The disk wind has a velocity of 
Vwind ~ 500 km and is most likely launched at outer parts 
of the disk, rwind ~ 10'''"g (lUeda et al.l 12009), by irradiation 
from the central source. During the limit cycle, the absorp- 
tion lines are the deepest in the oscillation-low state and be- 
come weaker as the hard X-ray flux increases. This strongly 
supports the photoionization of the wind; the corresponding 
ionization parameters estimated from the population ratio be- 
tween Fe XXVI and Fe XXV are consistent with th e change 
of the luminosity (based on an XSTAR simulation; iKallmanI 
120031) . The self-shielding is not relevant for photo-ionization, 
as the wind we observe is located 20° above the disk plane. 
The wind is most probably developed by the strong X-ray irra- 
diation in the oscillation-low state, as indicated by its strong 
iron-K emission line. Once a disk wind is launched, it will 
travel over a time scale of rwind /^wind ^10"^ sec, and hence 
we can approximate that the same disk wind steadily exists 
during the limit cycles. 

In the stable state, we do not detect any iron-K absorption 
lines. This suggests three possibilities; (1) the wind is not de- 
veloped, (2) the wind is fully ionized, or (3) the scale height of 
the disk wind is too low to be observed. We consider the first 
possibility unlikely because the absorption fines fiave been de- 
tected in the same Class x state e ven when the ha rd X-ray flux 
is lower than in our observation ("Lee et al.l l2002h . Indeed, the 
presence of iron-K emission line suggests that illumination of 
the outer parts of the disk does occur, which can trigger the 
launch of the wind. The second possibility cannot be ruled 
out, if the density of the wind is much smaller than that devel- 
ops in the oscillation states. Then, the ionization parameter 
may become sufficiently large to be fully ionized, even if the 
luminosity is smaller than that in the oscillation-high state. 
The third possibility is also plausible. 

5.4. Origin of the Continuum 

We have shown that the broad band spectra in all the four 
states can be commonly explained by Comptonization of seed 
disk photons off optically-thick (r « 7- 10), low tempera- 
ture (Te w2-3 keV), non-thermal hyb r id plas mas. Similar re- 
sults are obtained by IZdziarski et al.l (120011) in another state 
of GRS 1915H-105 (Class 7 or State B). These findings in- 
dicate that the optically thick and low temperature Comp- 
tonization are common features in accretion flows at high 
Eddington fractio ns, as suggested bv iDone et al.l (120041) and 
lUeda et aH (12009 ^. In the stable and oscillation-high states, 
we find that the fractional power supplied to non-thermal elec- 
trons is quite large {Int/lh = 0.4-0.6) and the slope of the 
electron energy distribution is steep (Finj «3-4) compared 
with the results repo rted from other data of GRS 1915H-105 
(IZdziarski et al.ll2005i) . This may be related to the high lu- 
minosities in those states (2 x 10^^ erg s~'). Figure 7 plots 
the absorption-corrected EI(E) spectrum in the stable state 
with the separate contribution from the MCD, reflection, and 
purely "thermal" Comptonization component that would be 
observed without non-thermal electrons. This illustrates the 
importance of non-thermal processes in these data. 

With the eqpair modeling, we find that a direct MCD com- 
ponent is necessary in addition to the Comptonized emission 
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Fig. 7. — The EI(E) spectrum in the stable state corrected for the inter- 
stellar absorption and dust-scattering. The thick solid curve (red) represents 
the total, the dashed curve (cyan) the direct MCD component, the dot-dashed 
curve (black) the reflection component, and the thin solid curve (blue) ther- 
mal Comptonization component that would be observed without non-thermal 
electrons. 

in the stable and oscillation-high states, assuming the same 
temperature for the seed photons and MCD component. This 
implies that (Case I) the Comptonizing plasmas are patchy 
over the disk, or (Case II) they are located inside a "trun- 
cation radius" of the standard disk. In reality, such trunca- 
tion would be smooth; the disk may gradually change into the 
"disk-Hcorona" region, as suggeste d for the v ery high state of 
canonical black hole binaries (e.g.. lDone & K ubota 2006). In 
Case II, the assumption of the same spectrum for the direct 
MCD and seed photons do not hold, and a disk radius derived 
from such a simple model may be largely overestimated from 
the true disk size if the energetics coupling between the disk 
and corona is taken into account dPone & Kub ota 2006). 

From the observed flux and innermost temperature obtained 
from the fit, we can estimate the innermost radius of the 
accretion disk both for the direct MCD component (/?J^*^°) 
and that of the seed photons for the Comptonized component 
^^Comp^ We calculate /;^°™p by assuming the conservation 
of photon numbers in the Comptonization process, ac cord- 
ing to the formula (Al) of iKubota & Makishirnal (120041) with 
the left-hand term increased by a factor of two for the opti- 
cally thick case. The resultant values of R^^^ and /?|°'^' = 

y^(^MCD)2_,_(/jComp^2 ^re listed in Table 2, without any cor- 
rections for the boundary condition and color/effective tem- 
perature, which would increase the radius by a factor of « 1.2 
(IKubota et al.lll998h . In Case I, /;|°'^' corresponds to the true 
innermost radius corrected for the partial coverage by Comp- 
tonizing clouds. 

We find evidence for the evolution of ^|°'^' during the oscil- 
lation, which increases from Ri210 km in the oscillation-high 
state to ^400 km in the oscillation-med/low states, although 
the significance of the change is marginal due to the large un- 
certainties in the latter states. They correspond to k, lOrg and 
~ 20rg, respectively, by assuming the black hole mass of 14 
Mq. This implies that an inner disk forms as accumulated 
matter rapidly accretes onto the black hole in the flare phase. 
Such variation of the disk radius during limit-cycle oscilla- 
tions has been reported based on simple MCD + power law fit- 
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ting (e.g.. iBelloni etanil997bl: iMigliari & Bellonill2003h . al- 
though we obtain larger radii than these studies that refer only 
to the MCD component by ignoring Comptonized photons. 
In the stable state, the radius (/?|°'^' «420 km) is significantly 
larger than that in the oscillation-high state. Thus, if Case 1 is 
true, our results suggest that the standard disk does not extend 
down to the innermost stable circular orbit (ISCO) regardless 
of the black hole spin in all our states, unlike in some cases of 
the "soft state" of GRS 1915+105 (McClintock et al. 200^ 
iMiddleton et al.1 l2006t IVierdavanti et al.l |200 9). Such low- 
temperature/high-luminosity disk emission that leads to an 
apparently large disk radius is observed in th e very high state 
of other black hole binaries (iDone & Kubotall2006 ). when the 
Eddington ratio is higher than that in the high/soft state. At 
even higher Eddington fractions that trigger limit cycles, as 
in our case, the situation may be the same as in the canonical 
very high state. 

The Comptonization parameters also show significant 
change during the limit-cycle oscillation, besides the disk ra- 
dius (and disk temperature) as discussed above. In the os- 
cillation high state, the fractional power supplied to non- 
thermal electrons becomes the largest (/nt/^h ~ 0.4) and the 
electron temperature highest (Te = 2.8 keV), which decreases 
to Int/h ~ 0.2 and T^= 1.8 keV in the oscillation-low state, 
respectively. Thus, in the flare phase, very efficient particle 
acceleration must occur by converting the gravitational en- 
ergy of accreting mass. This may increase the temperature of 
thermal electrons and expand the size of the corona. As the 
corona becomes geometrically thick, the "self-shielding" ef- 
fect starts to work, preventing the Comptonized photons from 
irradiating the outer part of the accretion disk (i.e., little reflec- 
tion signals). A similar situation (geometrically thick corona) 
may be realized in the stable state as well, where the reflection 
strength is weaker compared with the oscillation low states. 

To summarize, with the simple spectra model presented 
above, we show that the dramatic change of X-ray flux and 
spectra during the limit-cycle oscillation can be explained 
mainly by (1) the evolution of the disk geometry and (2) that 
of physical parameters of Comptonizing corona, particularly 
the fractional power supplied to non-thermal electrons. We 
note, however, that the assumption of the MCD model (i.e., 
standard disk) in all states and that of the same temperature for 
the seed photons of Comptonization would be oversimplifica- 
tion. For instance, disk instability theories predic t that a slim 
disk appears in the oscillation high state ( Abram ^wicz et all 
[l988), which predicts different spectra from the MCD model. 
This is not ruled out from the current simple analysis. Nev- 
ertheless, we stress that the role of Comptonization becomes 
very important in determining the whole dynamics and en- 
ergy spectra of black hole accretion disks at high mass accre- 
tion rates, which must be taken into account in any theoretical 
studies. 



6. CONCLUSIONS 

1. We find that the Suzaku broad band spectra of GRS 
1915+105 band both in the stable (Class x) and oscil- 
lation states (Class 6) are commonly represented with a 
model consisting of a MCD model and its Comptoniza- 
tion with a reflection component, over which the opac- 
ity of the disk wind (including iron-K absorption lines) 
is applied. The Comptonized component dominates the 
flux above ~3 keV. 

2. The Comptonization is made by optic ally- thick (r « 
7-10), non-thermal / thermal {T « 2-3 keV) hybrid 
plasmas. The non-thermal electrons are produced with 
10-60% of the total energy input to the plasma with a 
power law index of ?a3^, which account for the hard 
X-ray tail above '--^50 keV. 

3. During the limit-cycle oscillation, the reflection 
strength, estimated by the iron-K emission line, is the 
largest during the dip phase but disappears in the flare 
phase. We interpret this as evidence for self-shielding 
effects that Comptonized photons are obscured by the 
surrounding cool region of the expanded disk when 
viewed at a very high inclination angle from the outer 
disk. The evolution of the disk geometry is in accor- 
dance with theoretical predictions. 

4. The disk wind, traced by iron-K absorption lines of 
Fe XXV and Fe XXVI, always exists during the oscil- 
lation and its ionization well correlates with the X-ray 
flux. This supports the photoionization origin. In the 
stable state, the iron-K absorption lines are not detected 
probably because it is highly ionized and/or the scale 
height is small. 

5. The disk parameters suggest that the inner disk struc- 
ture is similar to that in the very high state of black hole 
binaries. The spectral variability in the oscillation state 
is explained by the change of the disk geometry and of 
physical parameters of Comptonizing corona, particu- 
larly the fractional power supplied to the acceleration 
of non-thermal particles. 
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